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Abstract 

The general formalism of homogeneous nucleation theory is applied to study 
the hadronization pattern of the ultra-relativistic quark-gluon plasma (QGP) 
undergoing a first order phase transition. A coalescence model is proposed to 
describe the evolution dynamics of hadronic clusters produced in the nucle- 
ation process. The size distribution of the nucleated clusters is important for 
the description of the plasma conversion. The model is most sensitive to the 
initial conditions of the QGP thermalization, time evolution of the energy den- 
sity, and the interfacial energy of the plasma— hadronic matter interface. The 
rapidly expanding QGP is first supercooled by about AT = T — Tc = 4 — 6%. 
Then it reheats again up to the critical temperature T^- Finally it breaks 
up into hadronic clusters and small droplets of plasma. This fast dynamics 
occurs within the first 5 — 10 fm/c. The finite size effects and fluctuations near 
the critical temperature are studied. It is shown that a drop of longitudinally 
expanding QGP of the transverse radius below 4.5 fm can display a long-lived 
metastability. However, both in the rapid and in the delayed hadronization 
scenario, the bulk pion yield is emitted by sources as large as 3 — 4.5 fm. This 
may be detected experimentally both by a HBT interferometry signal and by 
the analysis of the rapidity distributions of particles in narrow p^-intervals at 
small \pt\ on an event-by-event basis. 

PACS numbers: 12.38.Mh, 24.10.Pa, 25.75.-q, 64.60. Qb 



Typeset using REVTeX 



* Alexander von Humboldt Foundation Fellow 



1 



I. INTRODUCTION 



The hadronization of quark— gluon plasma (QGP) possibly produced in the early Uni- 
verse or in ultra-relativistic heavy ion collisions has received much activity during the last 
decades [p]-[2Tll. Despite the significant progress in understanding of the variety of possible 
signals and features of the QGP, the nature of the phase transition (FT) between deconfined 
and confined phase is not clear yet. Assuming a first order PT, usually an adiabatic scenario 
is invoked to describe the conversion of plasma into hadrons. A few years ago, in [|I^ it 
was pointed out that the coarse-grained field theory of homogeneous nucleation p2|,^ can 



be relevant for the relativistic PT also. This scenario assumes the nucleation of hadronic 
bubbles, e.g. bubbles of pion gas, within the (initially homogeneous) supercooled metastable 
QGP as the starting point of the PT. These bubbles are nucleated because of the thermo- 
dynamic fiuctuations of the energy density in the system. Then bubbles with radii smaller 
than the critical radius Rc shrink, bubbles of the critical size are in metastable equilibrium, 
while bubbles with radii larger than gain in size and develop into the new phase. The 
treatment of the relaxation of the metastable state within the framework of the nucleation 
theory provides the fundamental nucleation rate [p3| -|25[], which expresses the number of 
viable nucleating clusters of the new phase via the equilibrium number of critically large 
clusters. 

Langer's theory has been applied [0,113] to calculate the hadronization of rapidly expand- 
ing baryon-free QGP, produced in heavy ion collisions at RHIC and LHC energies. Although 
this approach seems to be more realistic than the idealized adiabatic PT, several questions 
remain open: These calculations fl^ have shown that the latent heat released during the 
plasma conversion is not sufficient to prevent the strong (20 % — 35 %) supercooling of the 
system. Thus, the rapidly quenched system leaves the region of metastability and enters 
the highly unstable spinodal region. Here the theory of spinodal decomposition might be in 
order to describe the further evolution of the fiuctuations leading to the break-up of that 
system. 

Secondly, since the critical radius drops quickly when the temperature is lowered, the 
bulk creation of the hadronic phase should begin [|TB| when the bubble radii are r < 0.8 fm. 
Finally, in the non-scaling scenario the bubbles grow independently on the total expanding 
volume. Then the completion of the PT will be significantly delayed []T^. Even within the 
scaling scenario the time necessary for the completion of the PT varies from 50 to 90fm/c 
(depending strongly on the numerical value of the surface tension). 

Hence, from the above one may conclude that either 

— the homogeneous nucleation scheme is inappropriate to describe the hadronization of 
relativistic systems or 

— that some important features of the first order PT are still missing. The situation 
would change if it turned out that the amount of plasma converted into hadrons had been 
underestimated in earlier works. 

Recently, the calculation of the dynamical factor n governing the growth rate of sub- 
critical bubbles was reexamined [Q, and the size distribution of bubbles in configuration 
space has been used to estimate the supercooling of rapidly expanding QGP [^. The latter 
plays an important role in the hadronization of plasma produced in relativistic collisions. 
When the critical radius of hadronic bubbles drops due to the rapid fall of the temperature 
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in the expanding system, the subcritical bubbles transfer to the region of supercritical sizes 
for these new conditions in the system. These bubbles then stop to shrink and start to gain 
size, thus increasing the total volume of the hadronic phase. 

In the present paper we study the effect of the bubble size distribution on the dynamics 
of the plasma— hadrons phase transition. The paper is organized as follows: the model used 
to study the QGP hadronization is described in Sec. ||. Section |T| reviews the formalism 
of nucleation theory. The evaluation of both dynamical and statistical prefactors appearing 
in hanger's theory are discussed. A coalescence-type model for further evolution of the 
nucleated hadronic bubbles is proposed. The role of the initial conditions (as well as effects 
of variations of the model parameters, i.e. the value of the surface tension, the minimum size 
of the pionic bubbles, the non-scaling regime and dilution factor, the prefactors, etc., on the 
relaxation of the metastable QGP is studied in Sec. |IV|. Section |V| presents the investigation 
of the finite size effects, the creation of long-lived states with metastable QGP and hadronic 
bubbles, and temperature fluctuations in the system. Finally, the results are summarized in 
the Conclusions. 



II. THE MODEL 

We consider a QGP produced in collisions of two heavy ions at RHIC or LHC energies. 
It is assumed that the plasma is thermalized soon. A wide range of initial conditions is 
studied, see Sec.^. The expansion and cooling is ruled by relativistic hydrodynamics. 

When the plasma cools below the critical temperature T^, a first order phase transition 
is initiated by the appearance of hadronic bubbles. The bag model equation of state (EOS) 
for the QGP consisting of gluons and massless quarks reads 

P. = ^(l6 + |%)T^-i? = ST^-5 (1) 

with the number of flavors, n^, and the bag constant, B. The EOS of the relativistic pion 
gas is 

Pk = ^T' = a,T' . (2) 
By imposing the condition Pg = ph at T = T^. one may find the critical temperature 

To (3) 

For a two-flavored QGP with B^'^ = 235 MeV ^ this gives = 169 MeV. 

The total energy density of the mixed quark— hadron phase, e, can be treated in the 
capillary (thin wall) approximation as a linear combination of the energy density of the 
hadronic phase, Ch, in the fraction of the volume occupied by hadronic bubbles, h = Vh/Vtot, 
the energy density of the QGP, e^, in the rest of the volume = (1 — h)Vtot, and the energy 
density of quark— hadronic interface, es = aSh/Vtot- 



S; 



e = heh + {l- h)eg + . (4) 

vtot 
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Here a is the surface tension of the interface between the two phases. The last term in 



(I) is usually disregarded jTSHlTl though its contribution to the total energy density may 
be comparable with the other two: the ratio of the interface between the phases to the 
hadronic volume is not too small. For instance, assume that all bubbles are of the same 
radius, \R\ = a. Then, at a = O.lTj^ 
Therefore, the surface term in Eq.@ 
hadronic clusters. 



28] the surface energy density scales as eg ~ a et 



may be omitted only for sufficiently large (r > 4fm) 
Similarly to the energy density, the total pressure also consists of three terms 



P 



hph + {1 - h)pq + J2ciPL 



(5) 



where the pressure of the spherical surface of radius Ri is given by the Laplace formula, 
p^"* = 2a /Ri, and Cj is the local concentration of bubbles of radius Ri in the total volume. 

The Bjorken model of scaling longitudinal expansion is applied to find the time 
evolution of the energy density. It yields the derivative of the energy density with respect 
to the proper time, r. 



de 
d7 



e + p 



T 



(6) 



Next one needs to compute the fraction of the total volume converted to hadronic phase. 
The pionic bubbles, appeared because of the fiuctuations in the energy density, will either 
shrink or grow. This problem cannot been solved on the basis of the thermodynamic theory 
of fiuctuations only. It requires a kinetic description of the process of bubble evolution. The 
volume fraction h{R, t) of bubbles of size R at time t obeys the equation of motion in the 
size space 



^^^ = -^HRMR,t)] + h-\R,t) 



(7) 



where v{R) is the radial velocity and h"''^'^'' denotes the hadronic fraction created at time t. 
Without the nucleation term Eq.(^) transforms into the continuity equation of the Lifshitz- 
Slyozov theory The volume fraction of bubbles of size R nucleated per unit 

time is given |3ll| by the distribution 



nucl 



(R) 



27r(9r + 2A|; 



: exp 



(8) 



X 



3r(a^ + a + 1) + 2A|a^ 



exp 



9r + 2A| 



1)^ 



da 



containing the nucleation rate /, the critical exponent r and the two dimensionless variables 
Xz and r (see below), so that 



POO 

/ K''"'\R)<1R = I . 
Jo 



(9) 



Instead of using a continuous spectrum, hadronic matter in our model is represented by a 
discrete spectrum of pion bubbles starting from tq = 1 fm. 
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III. THEORY OF HOMOGENEOUS NUCLEATION 



The homogeneous nucleation has been subject of intensive investigation both theoreti- 
cally and experimentally for a long period (for reviews see [p^-|3^ and references therein). 



In our analysis we follow the coarse-grained theory by Langer [^,^], who extended the 
classical Becker-Doring-Zeldovich (BDZ) theory of nucleation [^,0] to field theories. The 
nucleation rate in both classical and modern coarse-grained field theory reads 

/ = /oexp(-^). (10) 

Here Iq is the preexponential factor and AFc is the excess free energy of the critical cluster in 
the system. In hanger's theory the prefactor Iq is a product of a dynamical and a statistical 
prefactor, k and Qq, respectively 

lo = T^^o • (11) 
Zn 

It is interesting that under certain assumptions the prefactor derived in the classical theory 
may be obtained [^] identical to that of the modern theory. 



To clarify the meaning of both the dynamical and the statistical prefactors, let us consider 
a classical system with degrees of freedom described by a set of collective coordinates 
7], i = 1,. . . ,N. The coarse-grained free energy functional F{r]} of the system has local 
minima F{r]i} in the {r^j-space, corresponding to the metastable and stable states, separated 
by the energy barrier. The point of minimal energy along the barrier is the so-called saddle 
point {r]^}. Note that this saddle-point configuration corresponds to the critical cluster of 
a condensing phase in the classical theory. In contrast, in field theory the critical cluster 
of a condensing phase may not necessarily be a physical object but rather corresponds to 
a certain saddle-point configuration in phase space. The phase transition occurs when the 
configuration {rji} moves from the vicinity of a metastable minimum to the vicinity of a 
stable one. When the potential barrier is overcome, it is most likely for the trajectory of 
the system to pass across a small area around the intermediate saddle point {rj^}. The rate 
of the decay of the metastable state is determined by the steady-state current across the 
saddle point from the metastable to the stable minimum of F{ri}. 

Performing the Taylor series expansion of F{t]} around the {t]^} and keeping only 
quadratic terms yields 

1 ^ 

F{r^}-F{v'} = -Y^X,{v,-vl)\ (12) 

fe=i 

where are the eigenvalues of the matrix Mij = d"^ F / driidrij , evaluated at the saddle 
point. By definition, at {rj^} the free energy density functional reaches its local maximum. 
Therefore, at least one of the eigenvalues Xk must be negative. Following Langer, we denote 
it as Ai. Then one may approximate the potential barrier between the metastable and stable 
well by the excess of the Helmholtz free energy, AF, corresponding to the formation of a 
spherical bubble of size R. In the thin wall approximation, AF is the sum of the bulk and 
the surface energies. 
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AF{R) = -—R^Ap + inR^a , (13) 

with Ap being the difference in pressures inside and outside of the bubble. In the droplet 
model of Fisher the activation free energy includes also the so-called curvature term. 
It arises due to the small fluctuations in the shape of the bubble, which leave unchanged 
both, the volume and the surface area of the bubble, 

4 R 

AF^{R) = — TrR^Ap + AnaR^ + 3rT In — . (14) 
3 ro 

Here r ^ 2.2 is the Fisher critical exponent and tq is the radius of the smallest bubble in 
the system. Minimization of AF with respect to the radius R yields the free energy of the 
critical bubble 

AFf = lvri?,V + rT(31n^-l) . (15) 
Here the critical radius Rc should be evaluated by solving the equation dAF^ /dR = 0. 



It is convenient to introduce new variables [^]: the similarity number, Xz = Rc y^^ira/T, 
and the reduced radius, r = R/Rc. In terms of these variables we have 

AFF^ a 7 „ , R^ ,^ „v 

_^ = -, + ^ + 3rlu-^ (16) 



^ r+-Ai r^ + Air^ + 3rln^. (17) 

T V 3 / ro 

In the harmonic approximation for the activation energy of a bubble near the critical 
radius, Eq.(p!7|) reads 

AF^ 9r + 2\\ 



T 2 

and we get finally for the only negative eigenvalue Ai: 

Ai = -T (9r + 2A|) . (19) 

This expression will be used to determine the nucleation rate of the process, in particular 
the statistical prefactor. 



A. Dynamical prefactor 



The dynamical prefactor, k = d/dt [ln(i? — Rc)]-, which is related to the single negative 
eigenvalue Ai, determines the growth rate of the critical bubble of size Rc at the saddle 
point. To compute k one has to solve the hydrodynamic equations p^,|40[| which describe 



the growth of a bubble of the hadronic phase due to the diffusion flux through the interface. 
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For a baryon-free plasma, where the thermal conductivity is absent because of the absence 
of a rest frame defined by the baryon net charge, the dynamical prefactor has been calculated 
by Csernai and Kapusta [O to 



containing the bulk and the shear viscosities, ( and rj, and the difference Auj between the 
enthalpy densities of the plasma and hadronic phase, u = e + p. 

Here it is implied that the energy flow, uv, where v is the velocity of the net particles, 
is provided by viscous effects. Recently Ruggeri and Friedman p6[ argued that the energy 
flow does not vanish even in the absence of heat conduction. Since the change of the energy 
density e in time is given in the low velocity limit by the conservation equation, 

1^ = -V ■ M , (21) 

this means that the energy flow oc ujv is always present. Then the calculation of the dynam- 
ical prefactor for a system with zero thermal conductivity leads to the expression 

which agrees with the result of |^ obtained for nonviscous systems, i.e. the viscous effects 
cause only small perturbations to Eq. p^). Note that the dynamical prefactor K2 might 



violate the dynamical scaling laws of Kawasaki |Q in the vicinity of the critical point (for 
details see |]iO|Ji^ , |i3[] ) and, therefore, should be handled with great care. 

Figure |I] shows the evolution of the dynamical prefactors, ki and K2, below the critical 
temperature for different values of the surface tension. To calculate ki we use the fact that 
the shear viscosity of a two-flavored QGP, 77 = 1.29T'^/[a| ln(l/a5)] |^|, is much larger 



as compared with the bulk viscosity (. One can see that Eq. ( P2[ ) predicts higher rates n 
for moderate or weak supercooling. It means, in particular, that the minimal temperature, 
reached by the system during the cooling stage, should be higher if K2 is used in the cal- 
culations rather than ki. We will discuss the effect of the replacement of the dynamical 
prefactor on the course of the plasma hadronization in the model in Sec. fVj 

Homogeneous nucleation theory permits us also to determine the macroscopic radial 
velocity of bubbles via the dynamical prefactor, k, and the critical radius. Re- Assuming 



that bubbles of hadronic phase grow due to the diffusion flux through their interface we have 

viR)^^ = \K\ (^)\r-R^). (23) 



dt ' ' \ R 

U R = Rc, the radial velocity drops to zero and the bubble is in (metastable) equilibrium. 
If -R < Rc, v{R) is negative and, hence, the bubble collapses. If i? > Rc, the radial velocity 
is positive and the bubble grows. 
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B. Statistical prefactor 



The statistical prefactor, f2o, is a measure of the volume of the saddle-point region 
in phase space available for nucleation. Sometimes Qq is called also generalization of the 
Zeldovich factor Z although this is a crude simplification, since the relation between 
these two factors is actually more complex [R^. The product of Qq and the exponential 



exp (— AFc/T) gives the probability of finding the system at the saddle-point— rather than 
at the metastable configuration. 



According to [22-25], the statistical prefactor can be written as 




det (Mo/27rT) 
det (M'/27rT) 



1/2 



(24) 



where V is the available phase space volume at the saddle point, the index "0" denotes the 
metastable state, and the prime indicates that the negative eigenvalue Ai, as well as the zero 
eigenvalues of the matrix Mjj, is omitted. 

The calculation of the fluctuation determinant in Eq. (p^ ) is usually extremely difficult 
and, moreover, a very important uncertainty exists in the determination of Qq. Indeed, in 
the harmonic approximation (|12|) for the free energy density functional F{'r]} 



fin = V 



27rT 



1/2 



|Aii; 



( 



27rT 



N 

l=l0+2 \A; / 




(25) 



Here A['^^ and Xf are eigenvalues of the mobility matrix M , evaluated at the saddle point 
and at the metastable point, respectively, and Iq is the total number of symmetries of 
{F} which are broken by the presence of the saddle-point configuration. Since it is the 
translational symmetry of the system that is broken due to the bubble creation, the three 
translation invariance zero modes omitted in the fiuctuation determinant give rise to the 
prefactor proportional to A^ in the expression for the available phase space volume V: 



,(0) 



V = V 



3|Ai| 



3/2 



(26) 



Here V is the total volume of the system. In the fiuctuation corrections in the products 
over A^'^'' and A^^^* from ( ^5]) are absorbed into the free energy of the metastable region and 
the saddle-point region, namely 



exp i-Fo/T) ^ exp {-F{v'}/t) J] - 

1=1 \ \ 

exp(-F5/T)^exp(-F{ry^}/T) J] 



27rT 



1/2 



(27) 
(28) 



Therefore, the activation energy of a critical cluster is simply AF^ = Fs — F^ and 



fin = 1^ 



\i? vsj 



(29) 
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In Ref. [42 1 it was mentioned that there are four more terms in the product over Aj°^ than 



i s) 

in that over . Therefore, the free energy difference can not be precisely a logarithm of 
these products. The final expression for the statistical prefactor, which accounts for the four 
unpaired A^^^^'s modes, reads 




^o = -^^P^ ^ , (30) 



where ^ is the correlation length. 

The statistical prefactors corresponding to the expressions ( ^Of) and (|30D are shown in 
Fig. 0. In contrast to the curves in the lower panel, the curves in the upper panel diverge at 
T = Tc. This is an unphysical result, which is due to the fact that the nucleation process is 
turned off as T — > T^.It is worth noting, however, that the rise of Qq is counterbalanced by 
the factor exp (— AF^/T) that drops to zero at the critical temperature. Therefore, the total 
nucleation rate will increase with dropping temperature for both prefactors. Also, the cor- 
relation length, ^, is not a constant, but rather scales with the temperature in the proximity 
of the critical point, where the PT turns to second order, as ^{T) = ^(0) (1 — Tcrit/T)~" , 



with the critical exponent v = 0.63 [^]. Since such a critical point does not exist in the 
plasma to hadron gas phase transition scheme presented here, Eq. (^91) is used to calculate 
the statistical prefactor VLq. 



IV. RELAXATION OF THE METASTABLE QGP 

Like in any phenomenological model, the scenario of a homogeneously nucleating QGP 
is based on a set of model parameters. Here, the role of each parameter for the dynamics 
of the PT is studied in order to reveal those most relevant parameters to which the model 
results are most sensitive. 

The theory of homogeneous nucleation is valid for systems which are not too far from 
equilibrium. In particular, the supercooling in the system should not be too strong — 
otherwise the nucleation theory fails! The minimum temperature reached by the expanding 
and cooling plasma is shown in Fig. |^ as a function of the surface tension, a, and the critical 
time, Tc- The most reliable value of the surface tension lies within the range 0.015 < a /T^ < 
0.1 i.e. 2 < a < 12.5MeV/fm2 for the given T^. Note that even in the case of an 

extremely fast expansion and very low values of a, the supercooling of the system does not 
exceed 8%. This result does practically not change when bubbles grow independently in the 
non-scaling regime on the growing total volume [Fig. 0(b)]. Even when the surface tension 
and pressure of the bubble surface is included in the expressions for total energy density and 
pressure, Eqs. (§) — (|), [Fig. |5(c)], the results stay put. 

If the dynamical prefactor given by Eq. (|20|) is replaced by that of Eq. (|2^), the amount of 
the QGP volume converted into hadrons at the earlier stages of the supercooling is increased. 
Thus, the supercooling of the system in the latter case is about 2% weaker [Fig. |5(d)] than 
that shown in Fig. |^(a) — (c). 

In the analysis of the QGP equilibration times within the parton cascade model (PCM) 
by Geiger [^, which has been done for the RHIC and the LHC energies, yields the earliest 
equilibration time for the plasma as Tinit = 3/8fm/c with Tjnit ~ 2Tc. Then the critical 
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temperature will be reached at Tc = 3fm/c, which corresponds to a supercooling of about 
5—6% in our figures. Therefore, one may conclude that the application of the homogeneous 
nucleation theory to the hadronization process of even a relativistic expanding QGP seems 
quite reasonable. 

The temperature is plotted as a function of the proper time r in Fig. ^ for Tc = 
1.5, 3 and 6fm/c. The upper panels in this figure corresponds to conditions (a) of Fig. |^, 
the middle row corresponds to conditions (b), and the lower row corresponds to conditions 
(c) of the same figure. Hadronization causes the release of the latent heat, and the system 
reheats to temperatures close to T^. Then the nucleation and growth of hadronic bubbles 
comes to a halt. 

The continuing increase of the total volume leads again to a decrease of the temperature, 
and the phase transition continues immediately. These oscillations of the temperature in 
the vicinity of are well observable in Fig. ^. The mixed system is quite unstable at 
this stage, since any negligibly small rise of the temperature forces the system to reach the 
critical point, where, in turn, it may break-up into fragments. Then the theory of spinodal 
decomposition might be relevant to describe the hadronization of the rest of the QGP. The 
time At needed for the system (for cr = 5 — 10 MeV/fm^) to reach the zone of the oscillations 
scales with Tc as 

Ar = 2.9r,^/% (31) 

within the interval 1.5 < < 9fm/c. 

Changing the dynamical prefactor from ki to k,2 (Fig-i|) leads to negligibly small shifts in 
the time needed to reach the vicinity of the critical temperature. This is due to the fact that 
the late time evolution of the system is governed by the Lifshitz-Slyozov dynamics, which is 
generally much more important [^] for the course of a first order PT than the initial size 



distribution of clusters given by the nucleation theory. 

Therefore, the model does not appear to be very sensitive to a non-scaling growth of 
the nucleated bubbles, to the incorporation of the surface entropy in the rate equations, 
and to the numerical values of the dynamical prefactors obtained within the range of model 
parameters applied. The most important parameters are the value of the surface tension a 
and the time to reach of the transition temperature Tc, which is determined by the initial 
conditions and by the expansion dynamics of the system. The effects of varying just these 
two factors for the relaxation process of the metastable QGP are studied below. 

From Fig. |^ one may conclude that, at the very beginning of the phase transition, the 
process of bubble nucleation is the main mechanism of plasma conversion. The creation of 
the new phase reaches its saturation value soon after 0.5— 1.5fm/c, and the growth of the 
total hadronic fraction of the total volume proceeds due to the diffusion growth of already 
nucleated bubbles. Both the nucleation and the diffusion process contribute (almost equally) 
to the hadronization of the QGP. 

The critical radius in the system initially drops to about Ifm (Fig.|^), then, due to 
reheating, it rises up to i?c ~ 2 — 4 fm. Then the oscillations begin. At T = Tc, the critical 
radius is, of course, infinite. The weak changes of the temperature near the critical point 
cause significant oscillations in the value of the critical radius which, in turn, are responsible 
for the irregularities in the size distribution of bubbles at the final stage of the PT (see Fig. |^, 
right lower panel). The radial oscillations of the intermediate-sized hadronic bubbles can 
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result in the pulsed emission of matter and radiation owing to the strong acceleration and 
deceleration of matter in the vicinity of the bubble surface. These oscillations are analogous 
to the pulsations of a hot quark blob, discussed in |2^. Note also, that the bubble pulsations. 



occurring in medium, generate the sonic waves in the expanding plasma, but the discussion 
of this topic lies out of scope of the present paper. 

When the nucleation just begins, the size distribution of hadronic bubbles shown in Fig. |^ 
has a characteristic plateau-like profile. At Ar = Ifm/c, just after entering the metastable 
region, the critical radius drops to its minimal value of about 1 fm. Then practically all 
bubbles are growing. Reheating of the system leads to the rise of the value of the critical 
radius. As a result, a noticeable fraction of shrinking bubbles appears. The central plateau 
becomes narrower, also because of the increase of bubble density per unit of radial interval. 
At the end of the homogeneous nucleation stage, as mentioned above, the dips and peaks 
in the bubble size distribution arise due to the temperature oscillations. Note that the 
bubble size distribution established in the first order PT deviates clearly from the power- 
law distribution f{A) oc A"'^ {A being the size of a cluster), which is typical for the second 
order PT. 

The scaling of the change of the average radius of the hadronic bubbles with time (Fig. ^ 
demonstrates that the average radius depends mainly on the duration of the phase transition, 
but not on the expansion scenario. Since the nucleation of new bubbles is turned off after 
To ~ 1.2fm/c, we fit the distribution to the power-law 

< i?(r) > - < R{to) > = const x (r - tq)^/^ (32) 

with R{tq) = 3 fm, which is the Lifshitz-Slyozov (LS) t^/'^— law pO] of the coalescence process. 



We see that at the coalescence stage of the plasma conversion the agreement with the LS 
law is good. The deviations from this law at the late stage of the hadronization of QGP are 
also caused by the temperature fluctuations. 

Finally, the dependence of the results on the minimum radius of the nucleated bubbles is 
studied. The upper panel of Fig. |10| depicts the temperature curves calculated with vq = 1 
and 2fm, and the lower panel compares the size distributions of the bubbles at the end of 
the PT. Again one can see that the effect of the cut-off of small bubbles is negligible because 
of the rather broad initial distribution of bubbles in size space. 

The results of this Section may be summarized by concluding that the most important 
parameters of the model are the initial conditions, time, Tinu, and temperature, Tinu, of 
the QGP thermalization, as well as the scenario of further plasma expansion, and the value 
of the surface tension a. All other factors cause only small deviations from the solutions 
obtained with the fixed set of the aforementioned parameters. Note that effects like the 
final transients near the critical point are sensitively dependent on the fluctuations in the 
system, and these might (or might not) wash out these transients. Real systems produced in 
experiments with heavy ions are not infinite. Therefore, to complete our analysis, we have 
to clarify the role of the finite size effects. 



V. ROLE OF THE FINITE SIZE EFFECTS 

The biggest problem which has to be overcome before a careful investigation of the finite 
size effects can be tackled is the problem of the first order phase transition itself. Then, a 
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consistent treatment of the finiteness of the system and, especially, of the surface induced 
phenomena is still an open question. Attempts to estimate the role of the finite size effects 
on the phase transitions in nuclear matter have been made in PUI-HM] on a basis of a purely 



thermodynamical picture of fluctuations. Can this approach be modified and applied for our 
kinetic analysis of the plasma hadronization? To answer the question note that the situation 
we face here may be subdivided into two cases. 

First, let us imagine that a droplet of plasma is immersed into a hot gas of hadrons which 
acts as a heat bath, i.e. the temperature fluctuations in the droplet are suppressed. The 
compound system expands longitudinally, and the final size effects come into play via the 
finiteness of the transverse direction. The results obtained in Sec. |V| are valid for central 
collisions of gold or lead ions {R > 7fm) at relativistic energies. If then, by chance, the 
transverse radius of the expanding cylinder will be smaller than 7fm, the cut-off of large 
bubbles should reduce the volume fraction occupied by hadrons and affect the course of 
the phase transition. Figure ^ presents the evolution of both the temperature and the 
bubble size distribution as calculated with maximum radius 5fm and 4fm, respectively. 
The behavior of the system differs drastically between these rather close values of R: as the 
central plateau in the size distribution lies within the range of 3 < r < 4.5 fm at t ~ 10 fm/c, 
the cut-off of bubbles with i? > 5 fm does not cause noticeable deviations from the scenario 
discussed in Sec. However, for i? < 4fm the maximum temperature reached by the 
system during the reheating is reduced and the conversion of the QGP into hadrons is 
significantly delayed. As a result, the bubble size distribution has a pronounced peak at 
i? = 4fm (Fig. |TT|, lower panel). 

The values of the cut-off radius, at which the slowly varying temperature of the long- 
lived object does not exceed the 0.98 Tc upper limit, are listed in Table |for various critical 
temperatures, surface tensions and expansion rates. It is easy to see that all values are within 
the range of the central plateau in the bubble size distribution. Thus the most probable size 
of the emitting sources is not affected by the finiteness of the system. 

For radii smaller than 2 — 3fm, the surface induced effects must be taken into account. 
Undoubtedly, the theoretical treatment of the interface between plasma and hadronic matter 
is oversimplified as compared with realistic system. Therefore we are not able to make any 
quantitative predictions, based on the homogeneous nucleation theory, for the hadronization 
of such plasma filament. 

Second and perhaps more realistic is the case of the formation of a plasma pattern as an 
isolated system, i.e. without any contact with a heat reservoir. Although the energy of the 



expanding system is conserved, fluctuations of the temperature [4S] of the order of 



< i^Tf >= ^ (33) 

should occur. Here Cy = {dE/dT)y is the heat capacity at constant volume. For the two- 
component system consisting of QGP and hadronic matter, the temperature T is smeared 
around its mean value with the width 



The smaller the volume V of the system, the larger the temperature fluctuations, and vice 



versa. Moreover, from Fig. |T2| one may conclude that these fluctuations (at any given T) are 
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larger in the hadronic system rather than that in the pure QGP phase. This is due to the 
smaller specific heat of the hadronic phase. We have seen in Sec. |^ that — at temperatures 
of about 0.98 Tc (or lower) — the critical radius and nucleation rate depend only very weakly 
on temperature variations within the range AT/T^ ~ 2 %. 

The situation changes completely when the temperature approaches the critical one. 
Since the hadronic matter here occupies already about 70 % of the total volume, fluctuations 
of the temperature should be as large as 1 % of Tc, i.e. 1.7 MeV for V = 10'^ fm^. For smaller 
volumes, say 10 fm^, the fluctuations rise to 10 %, while for a volume as large as 16 ■ 10^ fm^ 
the width of the temperature smearing drops to 0.25 % of T^, still almost fully covering the 
range of the temperature fluctuations. Thus the expanding mixed system should break- 
up into fragments even earlier than in the scenario with the presence of heat bath. The 
conclusions drawn in Sec. |IV| remain true except for the appearance of irregularities in the 
size distribution of small and intermediate-size hadronic bubbles. 

VI. DISCUSSION AND CONCLUSIONS 

It has been shown that the theory of homogeneous nucleation is applicable to describe 
the hadronization of a relativistically expanding quark-gluon plasma produced in heavy-ion 
collisions. We have proposed a coalescence-type model to follow further the evolution of 
hadronic bubbles produced in a metastable QGP. The change of the average radius of the 
bubbles with time is shown to be consistent with the LS power law < R >(x t^^^. 

Various sets of model dependent parameters are used to study the role of each of them 
on the course of the plasma— hadrons FT. With rather good accuracy the number of these 
parameters may be reduced to the few main ones, namely initial conditions of plasma ther- 
malization, expansion scenario, and the value of the surface tension of the interface between 
plasma and hadronic matter. The supercooling of expanding QGF is found to be relatively 
moderate, 5 — 6 % only. Then the system reheats up to the critical temperature, where the 
temperature oscillations may occur. At At = 5 — lOfm/c after the beginning of nucleation 
process (this time depends strongly on the initial conditions and model of expansion) the 
system hits the critical point. Since at that time already about 70 — 80 % of the total volume 
is occupied by hadronic matter and the system is dilute compare to initial state, the rest 
of the QGF may not be sufficient enough to "glue" the hadronic bubbles in a compound 
state. The expanding system simply breaks up into fragments: hadronic clusters and small 
droplets of plasma. Due to the finite size effects the hadronization of QGF may be delayed 
and the long-lived objects containing plasma and hadronic bubbles are produced. But the 
temperature fluctuations cause the broadening of the critical temperature region and earlier 
disintegration of the system. One has to appreciate, however, that using the given EOS it 
is impossible to get the conversion of total amount of QGF faster than that of the idealized 
adiabatic transition. It means that the small sources like the QGF droplets will burst, emit- 
ting hadrons, from time to time up to about 40 — 50fm/c, while the bulk amount of plasma 
is converted into hadrons within flrst lOfm/c's. 

The entropy increase during the FT stage in the proposed scheme is small, 2 — 7%, and 
there should be no significant changes in the integrated particle yields or pion to baryon 
ratios which can be detected experimentally. On the other hand, the formation of large 
clusters and their further disassembly will lead to the significant multiplicity fiuctuations 



13 



in rapidity spectra of secondary particles which should also have low transverse momenta. 
To search for these clusters experimentally one can examine the rapidity and azimuthal 
distribution in small p^-intervals, looking for islands in sea of empty bins |]50[ . 

The nucleation process enforces the softening of EOS and diminishes the transverse flow. 
Then, the cluster size distribution in the model discussed departs from the simple power- 
law falloff, which is expected for the second order PT. The most probable radius of emitting 
sources, which is fully determined by the evolution of the value of critical radius with the 
temperature and by kinetics of the PT, varies from 3fm to 4.5 fm. This signal may be 
checked by the analysis of data on Hanbury-Brown-Twiss (HBT) correlations. However, the 
rest of the plasma dispersed between the hadronic bubbles is hadronizing also, giving rise 
to a substantial yield of small bubbles. Therefore, the presence of a plateau in the range 
i? = 3 — 4.5fm in the size distribution of hadronic clusters can be considered as a signal of 
the first order phase transition. 
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FIG. 2. (a): Statistical prefactor i7o given by Eq.(|30D vs temperature with (7 = 5 (solid line), 
10 (dashed line), and 15 (dotted line) MeV/fm^. (b): The same as (a) but for Eq.(|2|). 
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FIG. 3. Supercooling of the system as function of critical time Tc and surface tension a. (a): 
bubbles grow together with the total volume (scaling regime), entropy of the bubble surface is 
neglected, (b): bubbles grow independently on the total volume (non-scaling regime), entropy of 
the bubble surface is neglected, (c): non-scaling regime, surface entropy is taken into account, 
(d): the same as (c), but prefactor k\ is substituted by K2- 
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FIG. 4. The temperature evolution with time for the expansion scenarios with Tc = 1.5 fm/c 
[(a)-(c)], 3 fm/c [(d)-(f)], and 6 fm/c [(g)-(i)]. The sequence of conditions within each subgroup 
of panels corresponds to that of Fig.^ (a)-(c). The dotted curves are the idealized adiabatic 
scenario of the PT, the others show the result of solving of rate equations with a = 5 (solid lines) 
and 10 (dashed lines) MeV/fm^. 
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FIG. 5. (a) - (c): the same as Fig.|| (a),(d), (g) with a = 5MeV/fm^ but for the solutions 
of the rate equations with ki (sohd hnes) and with K2 (dashed hues). 
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FIG. 6. The part of the QGP volume converted into hadrons for the ideahzed adiabatic PT 
(dotted Unes) and for homogeneous nucleation scenario with cr = 5 MeV/fm^. SoUd lines correspond 
to the total volume fractions of hadronic matter, dashed lines denote the increase of hadronic 
volume due to the nucleation of new bubbles, and dash-dotted lines indicate the enlargement of 
the hadronic bubbles due to diffusion. 
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FIG. 7. Change of the critical radius as a function of time r — Tc for the expanding system 
with Tc = 1.5 (a), 3 (b), and 6 (c) fm/c, with cr = 5 (soUd hnes) and 10 (dashed hnes) MeV/fm^. 



23 





lU 




1 u 




10 




10 








10 


Cm 






10 




10 




10 




10 




10 



-4 
■5 
■6 
■7 
■8 
■4 
-5 
■6 
-7 
-8 



At = 1 fm/c 



■m- 



Ax = 3 fm/c 





2 4 2 4 

R (fm) 

FIG. 8. Size distribution of hadronic bubbles at Ar = 1,2, and 3 fm/c after the beginning of 
the nucleation, and at the freeze-out for r = 1.5, At = 4fm/c (solid curves), r = 3, At = 6fm/c 
(dashed curves), and t = 6, At = 11.5 fm/c (dotted curves). 
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FIG. 9. Average radius of hadronic bubbles as a function of r — Tc for the longitudinal 
expansion with a = 5MeV/fm^ and with Tc = 1.5 (triangles), 3 (circles), and 6 (squares) fm/c. 
Solid curve is the fit to Lifshitz-Slyozov power-law < R>oc t^l^. 
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FIG. 10. (a): Temperature as a function of time for the expansion with a = 5MeV/fm^ and 
Tc = 3 fm/c. Curves show the result of solving of rate equations with minimum radius of hadronic 
bubbles Rmin = 1 (solid line) and 2 (dashed line) fm. 
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FIG. 11. Upper : The same as Fig. 10 but for maximum radius of hadronic bubbles -R^ 



(solid line) and 4 (dashed line) fm. Lower : Size distributions of hadronic bubbles at the freeze-out 
for the SQ^niG scGna,rio cis in upper panel with Rrnax — 

5 (solid line) and 4 (dashed line) fm. 
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FIG. 12. Temperature fluctuations (< (AT)^ >)^/^/T as function of reduced temperature 
T/Tc and hadronic fraction h of the total volume of the system. 
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TABLES 



TABLE I. Maximum transverse radius of longitudinally expanding QGP at which the 
long-lived metastable state appears. The maximum temperature reached by the system during 
the plasma conversion is less than 0.98 Tc. 



Surface tension, 


Tc = 1.5 fm/c 


Tc = 3 fm/c 


Tc = 6 fm/c 


a (MeV/fm2) 


Rmax (fm) 


Rmax (fm) 


Rmax (fm) 


Tc = 150 MeV 








2 


3.64 


3.96 


4.35 


5 


3.70 


4.25 


4.63 


10 


3.85 


4.35 


4.76 


Tc = 170 MeV 








2 


3.40 


3.71 


4.05 
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3.60 


4.00 


4.35 


10 


3.76 


1.20 


i.50 


Tc = 200 MeV 








2 


3.13 


3.37 


3.69 
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3.23 


3.60 


3.98 


10 
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3.79 


4.20 



29 



